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NQM&NCLATORB 

A-  cross-sectional  area,  ft2 

Aw-  wall  area,  ft2 

C-  speed  of  sound,  ft/see. 

Cp*-  specific  heat  at  constant  pressure,  Btu/lbm  °R 

Cp-  specific  heat  at  constant  pressure,  ft  lbf/slug  °R 

D-  pipe  diameter,  ft. 

0-  mass  Telocity,  w»/A  lbm/sec  ft2 

g-  gravitational  acceleration,  32.2  ft/ sec2 

«c~  gravitational  constant  32.2  lbs  ft/lbf  sec2 

f-  friction  coefficient 

h-  enthalpy  per  unit  mass,  ft-lbf/slug 

J-  Mechanical  equivalent  heat,  77d  ft.lbf/B 

k-  ratio  of  specific  heats 

M-  Mach  Number,  v'JkRT 

P-  pressure,  lbf/ft2 

Q-  heat,  Btu 

B-  gas  constant,  Btu/lbm  °R 

R»-  gas  constant,  ft.lbf/lbm  °R 

T-  temperature,  °R 

w'«  mass  rate  of  flow,  lbm/sec. 

w-  mass  rate  of  flow,  slug/sec. 

V-  velocity,  ft/see. 

2-  elevation,  ft. 

Vg»-  velocity  component  of  gas  injection  in  X-direction,  ft/sec. 

v-  specific  volume,  ftVslug. 

Vg-  velocity  of  gas  injection,  ft/sec. 


u 


Greek  Letters 

p*  density,  lbm/ft^ 

p   density,  slug/ft^ 

7J/  well  sheering  stress,  lbf/ft2 
(o)  signifies  the  stagnation  state 
(  )*  signifies  state  at  which  the  Mach  Number  is  unity 

A   difference  between  values 


IKTRODUCTIOK 

The  usual  treatment  of  the  compressible  fluid  flow  is  to 
•elect  one  of  the  flow  parameters  as  the  independent  variable, 
and  to  determine  how  the  dependent  variables  change  when  the 
Independent  variable  is  allowed  to  vary.  The  usual  independent 
parameters  for  negligible  elevation  effects  are  area,  stagna- 
tion temperature,  friction  (which  causes  the  impulse  function 
to  change),  and  the  mass  of  the  fluid  flowing  through  the  duct. 
By  allowing  only  one  independent  parameter  to  be  changed,  dif- 
ferent kinds  of  the  compressible  flow  will  be  determined. 

In  the  case  of  simple  area  change  with  constant  stagnation 
temperature,  constant  mass  flow  through  the  duct  and  negligible 
friction  at  the  wall  of  the  duct,  the  flow  will  approximate  a 
reversible  adiabatic  process.  Here  it  is  said  the  flow  will 
approximate  a  reversible  adiabatic  process  because  in  actual 
practice,  if  the  variable-area  duct  is  long  and  poorly  insula- 
ted, the  process  will  not  be  adiabatic.  Where  there  is  a  poorly 
insulated  pipe  and  one  that  is  sufficiently  long  the  surf see 
area  will  be  large  enough  to  cause  heat  flow,  which  la  the  case 
for  natural-gas  pipelines.  In  the  case  of  a  well  insulated  long 
pipe,  or  an  uninsulated  short  pipe,  the  flow  will  approximate 
the  reversible  and  adiabatic  process.  Because  of  insulation 
there  will  be  negligible  heat  flow  in  a  long  pipe,  and  in  the 
case  of  an  uninsulated  short  pipe  the  area  for  heat  transfer 
will  be  small  and  the  heat  flow  will  be  of  a  negligible  amount. 
In  both  caees  the  pipe  la  assumed  to  be  perfectly  smooth.  The 


treatment  of  reversible,  adiabatic,  horizontal,  compressible 
fluid  flow  is  usually  found  in  text  books  on  gas  dynamics. 

Compressible  fluid  flow  in  ducts  having  constant  area,  in 
which  friction  is  present,  and  having  constant  stagnation  tem- 
perature and  constant  mass  flow  of  the  fluid  through  the  duct, 
will  give  rise  to  the  Fanno-Line  process.  In  this  type  of  flow 
there  Is  negligible  heat  transfer  and  the  system  approximates 
an  adiabatic  process.  In  the  case  of  a  long  natural-gas  pipe- 
line there  is  sufficient  area  to  cause  heat  transfer,  which  will 
give  rise  to  a  dlabatic  process.  For  short  ducts  the  area  for 
heat  transfer  is  small,  and  flow  of  heat  will  be  of  negligible 
amount.  For  the  actual  pipeline  which  is  well  insulated  the  flow 
will  approximate  the  Fanno-Line  process.  The  limiting  ease  for 
the  long  duct  with  frictional  flow  is  one  in  which  there  is  heat 
flow  and  there  is  constant  temperature  along  the  duct.  Actually 
in  the  Fanno-Line  process  the  one  independent  parameter  is  length, 
and  the  friction,  which  results  from  this  independent  parameter, 
causes  the  impulse  function  to  decrease.   This  type  of  flow  is 
discussed  in  most  of  the  text  books  on  gas  dynamics. 

The  case  of  diabetic,  constant  area,  constant  mass  flow, 
horizontal  flow  without  friction  is  termed  the  Rayleigh-Line  pro- 
cess. Such  ideal  conditions  do  not  exist  in  practice,  because 
the  frictional  effects  in  heat  exchangers  are  not  negligible. 
If  the  stagnation  temperature  is  raised  by  some  chemical  combus- 
tion process,  there  will  be  a  change  in  the  chemical  composition 
of  the  gas  flowing  through  the  duct.  If  this  situation  is  over- 
looked, and  if  the  flow  is  analysed  by  considering  the  simple 


stagnation  temperature  changes,  some  important  conclusions  for 
practical  significance  will  be  drawn.  This  type  of  flow  is 
called  the  Rsyleigh-Line  process,  and  it  is  treated  in  most  of 
the  books  on  gas  dynamics,  in  which  the  effects  of  elev  tion  are 
not  considered. 

In  this  report  analyses  were  made  on  flows  in  which  there  is 
no  heat  transfer  and  in  which  the  cross-sectional  ares,  the  im- 
pulse function  and  the  stagnation  temperature  are  all  kept  con- 
stant. In  the  first  section  an  analysis  was  made  in  order  to 
establish  the  relations  of  pressure,  density,  temperature  and 
mass  flow  rate  in  terms  of  Mach  Nunbor.  These  relations  are  in- 
dependent of  elevation  changes  and  friction.  From  these  formulas 
a  table  and  a  graph  of  the  properties  have  been  prepared  for  this 
fluid  flow. 

In  sn  actual  case  of  this  kind  of  flow,  the  frictional  ef- 
fects slong  the  wall  of  tht  duct  should  be  considered.  This 
actual  fluid  flow  was  treated  in  the  second  section  of  this  re- 
port, and  the  effects  of  Mach  Number  on  the  length  of  the  duct 
were  investigated  and  a  table  and  a  graph  have  been  prepared. 
In  this  investigation  the  gas  was  injected  with  its  forward  com- 
ponent of  velocity  equal  to  that  of  the  free  stream. 

The  final  investigation  was  msde  considering  the  effects  of 
elevation  and  neglecting  the  fluid  friction  at  the  wall  of  the 
duct.  First  an  analytical  investigation  was  made,  and  this  in- 
vestigation was  checked  with  numerical  calculations.  The  re- 
sults of  these  calculations  were  plotted  on  graphs. 

The  fundamental  equations  of  gas  dynamics  and  thermodynamics 


were  used  in  this  investigation.  The  differential  form  of  the 
equations  provides  the  information  required  to  show  the  signif- 
icance of  elevation  changes.  For  the  flow  having  constant  area, 
constant  Impulse  function,  and  constant  stagnation  temperature 
all  possible  states  are  realized  by  the  use  of  the  continuity 
equation,  energy  equation,  equation  of  state  and  momentum  equa- 
tion. 

SECTION  I 

Derivations  of  Pressure,  Density,  Temperature 
and  U»bb  Rate  Flow  in  Terms  of  Mach  Number  for 
Constant  Area,  Constant  Impulse  Function  and 
Constant  Stagnation  Temperature  Fluid  Flow 

For  this  analysis  the  stagnation  temperature  is  maintained 
constant  by  (1)  adiabatic  flow  with  the  stagnation  temperature 
of  the  injected  gas  equal  to  the  stagnation  temperature  of  the 
main  flow,  or  (2)  diabatic  flow  with  the  stagnation  temperature 
of  the  injected  gas  adjusted  to  maintain  a  constant  stagnation 
temperature  along  the  duct. 

The  Injected  gas  is  assumed  to  mix  completely  and  instan- 
taneously with  the  main  stream  flow,  i.e.,  the  mixing  distance 
along  the  duct  is  sero.  In  the  case  of  horizontal  flow  with  no 
friction  the  injected  gas,  in  order  to  maintain  a  constant  value 
of  the  impulse  function,  must  enter  the  region  of  the  main  flow 
perpendicularly.  In  the  case  of  elevation  changes  the  main  flow 
must  be  upward  as  will  be  shown  in  Section  III,  and  the  injected 
gas  must  have  a  component  of  velocity  in  the  upward  direction. 
In  the  case  of  horizontal  flow  with  friction  the  injected  gas 


oust  have  a  velocity  component  in  the  uirection  of  the  main 
stream  flow,  aa  that  will  be  shown  in  Section  II.  No  restric- 
tions are  pieced  on  the  transverse  velocity  of  the  injected  gas 
by  the  flow  of  heat.  However,  the  main  stream  flow  is  assumed 
to  be  one-dimensional,  i.e.,  any  cross- section  of  the  pipe  has 
uniform  properties,  and  the  mixing  length  along  the  pipe,  as 
stated  before,  will  be  assumed  zero. 

Impulse  function.  The  impulse  function  is  defined  ss 

F  z   PA  ♦  fV2A  «  - (1) 

for  a  perfect  gas,  and  from  the  equation  of  state, 
V2  -     P     v2     and  kRT  «     C2 

PV2  s     PkM2 (2) 

From  equations  (1)  and  (2) 

F«PA+  PkM2A  a  PA(l*kM2) 

T*  s  mfMlel***),     A  s  A*,     M*  =  1,  and  F*  s  F 

Temperature  ft*lflUoffilhlR-      From  the   steady  flow  energy  aqua- 
tion  and  the  requirement  of  constant  stagnation  temperatures, 

To  =  To     =  T<X  *  kd  M2)   s     T*(l  ♦  k^i  H»2) 

«  'frik+l),     H+  n  1 
2 

Density  Ratios.     From  the  equation  of  state 


From  equations  (3),    (4),  and  (5) 


4.  .,1       i*Y»?-2U*  ^"2) 


(6) 


s!Hfe»*    "    *fh      <7) 

From  aquations  (4)  and  (7) 

J-    =    r  /  ^  (d) 

2 

Ka;=s  ^ate  |  low    latio.      From   the  equation  of   continuity, 

fr   --flfer-        A   ■    A* 

From  aquations  (9),   (8),  and  (6) 


2 


"^  -  TuETPT     J  — Ei 

2 


JL     =     M>l2(k»l)    U»^M2)        (10) 

**  1  +  kM2 

In  ordar  to  determine  the  manner  in  which  the  Hach  Number 
may  change,  the  relation  between  entropy  change  and  Mach  Number 
change  will  be  investigated. 


Sfttropy  RfltttMiMp*     The  entropy  change  from  (    )*,  the 
state  where  K  =  l»o,  to  a  given  state  Is 


S*  -  S  *  CD     in 


-fc- 


.     R     In  Po* 


o  *o 

S  -  S*     -  -Bin     JP 

P* 
ro 


s^   . 


In 


From  the  isentropic  pressure  relationship, 


(11) 


-p2     =     <x  *  Jfel  M2) 


k/k-1 


2j    .     (l.J^2)     k/^ 
=     (1  *  lEii)         k/k-1» 


M*     «     1 


k/k-1 
Prom  equations  (12)  and  (13) 


*    Ujl) 


k/k-1 


*! 


Prom  equations  (14)   and  (3) 


'it-'  -H*» 

o 


ljk  M2) 


^¥«2 


-,  k/k-1 


k4l 


(12) 


(13) 


P.  P  (1  ♦  *gl  M»)     ^I1 (u) 

k/k-1 


(15) 


From  aquations  (11)  and  (15) 

S-S*  .  R       In  (l»kH2) 
(ilk) 


S-§£-     in     (jf*  M2   ) 
^T  (i4k) 


k+1 
2 


nkA-1 


1  ♦ 


^M2 


r    1** 


(16) 


All  tha  derived  relations  in  thia  aaction  ara  expressed  in  terns 

of  ttaoh  Number  only,  and  thua  ara  independent  of  elevation 

changes  and  friction. 

Differentiating  aquation  (16), 

k/k-l     r       ir    k4>1     _,k/k-l 

(iikiSi) 


i+ 


¥»2 


i*^ 


Irk 


^i*i=l  M2_ 


°!&* 


<<^) 


k/k-l 


¥  J 

1+k    M2 
(la  ^j**2  ) 


14k 


-,  k/k-l 


L^l 


k/k-l 


1+k  M' 


1+  ^ 


-,  k/k-l 


(1*  *jfe2)         "     (k)   d(K)2 


II±kK 


2)kIX   {*»  ^M2) 


iA-i 


f-sfcr  <,("'2-  !    • 


•  ¥ 


d  (K): 


1+ 


ifc 


TT2 


d  (K) 


^       2(1*11*)      (i^H2) 


(17) 


From  the  second  law  of  thermodynamics  mn  adiabatic  process 
can  never  be  accomplished  with  a  decrease  in  entropy.     An  ex- 
planation of  equation   (17)   shows  that 

(i)     In  the  case  of  M<1,   the  value  of  d(M)2  oust  be 
positive. 

(ii)  In  the  case  of  M>1,  the  value  of  d(M)2  must  be 
negative. 

It  is  concluded  that  for  subsonic  flow  the  Mach  Number  must 
increase,  anu  that  for  supersonic  flow  the  Mach  Number  must  de- 
crease. 

Formulas  For  Numerical  Calculation 

*L.  z      <VlAA  but  Ai.  A2 

w2  f2V2A2 

?2V2  *h  P2  ^5^2 

»    &     f    ft     e  /5 

P2  fc-2         V   TJ 

£*     «     ^     •     £*         (Th«  values  of  the  ratios  of  P/P*  csn  be 
f 2  P*  P2  taken  from  Rayleigh-Line  tables  because 

the  Impulse  function  is  constant  in  the 

Rayleigh-Line  process.) 

Tl     -    Ii     .     Io         (The  values  of  the  ratios  can  be  taken 
T%     "     T0  7~  from  the  isentroplc  tables.) 

Substituting  these  values  of  Ti   and  Pi   in  the  formulas  for 

v-2.  *     * 


"A    -  P 

w2      "   p* 


•«-7%'%-% 
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Table  1.  Properties  for  constant  area,  constant  impulse 
function  and  constant  stagnation  temperature 
fluid  flow  as  function  of  Mach  Number.* 


• 

Mach  t 
Number: 

• 
• 

P/P* 

• 
• 

:      T/T* 

• 
• 

:     <%* 

:     w/w*     : 

P°/Po*  J 

V/V* 

0.1 

2.3663 

1.1976 

1.9763 

0.21626 

1.2591 

0.1092 

0.2 

2.2727 

1.1904 

1.9090 

0.416597 

1.2345 

0.2160 

0.4 

1.9607 

1.1620 

1.6662 

0.72734 

1.1565 

0.4320 

0.6 

1.5957 

1.1194 

1.4255 

0.90494 

1.07525 

0.6346 

o.d 

1.2656 

I.O636 

1.1696 

0.96160 

1.0193 

0.6272 

1.0 

1.0000 

1.0000 

1.0000 

1.00000 

1.0000 

1.0000 

2.0 

0.36363 

0.66666 

0.54545 

0.69072 

1.5030 

1.6230 

3.0 

0.17647 

0.42657 

0.41176 

0.60666 

3.4244 

1.9620 

4.0 

0.10256 

0.26462 

0.35690 

0.76751 

6.2266 

2.1400 

6.0 

0.04669 

0.14634 

0.3 1906 

0.73234 

36.4950 

2.2960 

6.0 

0.02649 

0.06675 

0.30456 

0.71661 

136.6200 

2.3536 

10.0 

0.01702 

0.05714 

0.29767 

0.71205 

361.6100 

2.3900 

*Theee  values  are  independent  of  elevation  changes  and  frictional 
effects. 
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Fig.  1.  Properties  for  constant  area,  constant 

stagnation  temperature  and  constant  area 
fluid  flow. 
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SECTION  XI 

Analysis  of  Horizontal  Flow  With  Wall  Friction 

In  this  analysis,  the  affects  of  wall  friction  are  Inves- 
tigated through  the  introduction  of  the  friction  parameter. 
In  this  analysis  a  step-by-step  numerical  method  will  be  used, 
aa  no  exact  method  of  analysis  could  be  found.  The  steps  will 
be  divided  in  the  increment  of  0»1  Mach  Number. 
Assume:   1.  Flow  to  be  horizontal,  with  wall  friction. 
2.  One  dimensional  and  steady  flow. 

3*  Continuous  injection  of  gas,  with  complete  and 
instantaneous  mixing  of  the  injected  gas. 

4.  Constant  area. 

5.  Constant  impulse  function. 

6.  Constant  stagnation  temperature. 

7.  The  gas  is  perfect. 

8.  V_t  is  the  forward  component  of  the  velocity  Vg 

with  which  the  gas  is  injected.  Vg?  9   V,  the 
main  stream  velocity. 

Sample  calculations  to  illustrate  the  step-by-step  method 

of  analysis: 

*1  ■  0-5,  A  r  1  ft2,  Pi  r  100  psia,  Tx  =  1000  °R,  4f  s  «02  i.e., 

f  =  .005 

fl  =  £i    «  100  -     0-26995  lbm/ft3 

HTi         0.37043  x  1000 

vl  ■  Mlcl  s/kRTi'Mx  ,     k  =  1*4,     R  =     53*342 

Vl  s  0.5x49 «0172  /ioOO  =  775*08  fpa 
wi  s  Ai^Vi,  Ax  s  1  ft2 
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six  0.26995  x  775.08 
=     209.23  lbm/eec. 
M2     3     0.6 

T 

JL  s  0. 93284  (from  isentropic  tables),  T0  s  1050  °R 

To 

T2  s  1050  x  0.932*4  -  979.48  °R 

V2  s    M2C2    s  M2  /kRT2     «     0.6  x  49.017  x   y/ 979.48 

s     920.25  ft/eec. 
^.  1    1.5975     (Kayleigh  Line  tablet) 

£j_  =     1.7778,         P*     z     i  lift       s     56.24  p.U. 

P2  s     1*5975  x  56.24 
r    89.757  pel*. 

equation  of  state, 

^2     •    H         -    $9,759 - =       0.24739      Ba-'ft* 

STJ~  0.3705  x  979.48 

Continuity  Equation, 

*2     =     A2P2V2   ,     A2     s       1  ft2 

■  0.24739  x  920.25 
w2  s  227.66  lbm/see. 

w2"wl  =  18-23  lbm  -  ^wx 

Assume  the  length  of  the  pipe  for  thie  change  of  Mach 
Number  is 


14 


i  ! 


i — i 

(jy» ax *-@ 


Pig.  2 

X  fro*  section  (1)  to  section  (2)  then  the  wall  shearing  stress 
is  (Fig.  2) 

"Tw  -  f eyl dd) 

2 

In  the  present  analysis  an  averaging  method  will  be  used, 
therefore,  the  above  shearing  stress  will  be  reduced  to  sn  aver- 
age shearing  stress  between  section  (1)  and   section  (2),  and  it 
is  aasuaed  that  the  co-efficient  of  friction  is  constant  through- 
out the  flow. 

^     :     f^ 

T 

2  2x2 

«     0.26995  x   (775.08)2     ♦  0.21739(920. 25)2 

4x32.174 

*     2,  *90       (19) 

/  w       ■     -005  x  2090     (substituting  the  value  froa  (19)   in  (Id) 
=     14.45  lby/ft2 
The  wall  friction  between  Section   (1)   and  Section  (2). 

Fw     x     frictional  force,   -  shearing  stress  x  wall  area. 
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\,    a     wetted  area  of  the  wall. 
.    7T  jM.  A  -     1  ft2 

D     s     1.129  ft. 
A*     z    7T  x  1.12S&X. 

The  wall  frictional  force  is  set  equal  to  the  axial  compo- 
nent of  the  momentum  of  the  weight  of  the  gas  injected  to  get  the 
length  required  for  this  kind  of  change  of  flow  properties. 

4Fw=awV       (20) 

V  «  775.08  *   920.21     347.66  ft/sec. 

2 

4w    ■     lg,43       s     0.573  slug/eec. 
32.174 

F„     -     *-xl,129  x  14.45^ Xx     »     51.251  *i 

From  equation  (20) 

**!  x  51.25  «     0.573  x  847.66 

4^l-    ftt5734x  jilsH 
51.25 


s    9.465  ft. 


r3 


iSL    «    tostfltfsMl    =      166x10 
D  1.129 

These  calculations  were  repeated  for  other  increments  of  Mach 
Numbers.  The  results  are  plotted  on  Fig.  3,  in  which  the  "X"  in 
the  expression  ±fk   represents  the  distance  along  the  pipe  from  a 
particular  section  to  the  section    where  Mach  Number,  M  a  1.0. 

From  the  analysis  presented  in  Section  I,  the  flow  will  al- 
ways occur  in  such  a  manner  that  the  Mach  Number  tends  toward 
unity,  as  shown  by  the  arrows  on  the  graph. 


Fig.  3.  Dimensionlesa  pipe  length  parameter  for  horizontal     || 
flow  with  wall  friction*  kf&wx  vs  Mach  Number, 


£  ..&$::.io-.4\&0\Z--&:r1<A    i 


s^SBBffite  h-S  • 


iaatattaaaaaaama^^^^ttEc^^^iM ' ' ' ' ' ' ' ' ' ' ' im 


Table  2.  Dimenaionless  pipe  length  parameters  for 
horizontal  flow  with  wall  friction. 
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4fXg§ax  xio2 

1.0 

0 

0 

0 

1.1 

0.343 

0.343 

0.686 

1.2 

0.873 

1.216 

2.432 

1.3 

1.22 

2.436 

4.872 

1.4 

1.405 

3.041 

7.682 

1.5 

1.515 

5.356 

10.712 

1.7 

3.47 

8.026 

■ 

17.652 

2.0 

5.71 

14.536 

29.072 

1.0 

0 

0 

0 

0.9 

0.443 

0.448 

0.896 

0.0 

1.69 

2.138 

4.276 

0.7 

3.25 

5.488 

10.976 

0.6 

5.53 

11.018 

22.036 

0.5 

9.465 

20.483 

40.966 
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SECTION  III 


Mathematical  Analysis  of  Vertical  Flow  Without  Friction 


Vertically  Upward  Flow. 


1.  Constant  srea. 

2*  Constant  Impulse  function 

3.  Constant  stagnation  temperature. 

4*  Frietionless  flow. 

5.  Perfect  gas. 

6.  Continuous  gas  injection  with  complete  and  instan- 
taneous mixing  of  the  injected  gas. 

7*  One  dimensional,  steady  fluid  flow. 


a  t-l 


7 


SI 


/ 


© 


\ 


s 


\ 


Fig.  4. 

Consider  the  flow  in  a  duct  between  two  sectior s  as  in- 
finitesimal distance  dZ  apart.  In  this  element  of  duct  length, 
dZ,  gas  is  Injected  into  the  stream  at  mass  rate  of  flow  dw. 

Equation,  of  State. 

P  *  fRT 

iiE  =  4£  +  dl (21) 
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Definition  o£ 
IcST 

AiMii  ,  i*v£  .   di    mi 

Hi  v2  * 

Equation  of  Continuity.  The  injected  gea  flow  is  assumed  to 
be  injected  continuously  along  the  length  of  duct. 

It  is  assumed  that  the  injected  gas  and  the  main  stream  have 
mixed  perfectly  at  the  down  stream  boundary  of  the  control  volume. 
The  mass  flow  of  the  main  gas  stream  may  be  expressed  as  follows: 

w  s  PAY  ,    A  •  A*  «  constant. 

ifc  s  l£_  +  dV   (23) 

Energy  Equation.  The  energy  equation  for  the  flow  through 
the  control  surface  of  figure  may  be  written,  assuming  the  grav- 
ity effects, 

»4  j£_     *-m-  (2) 


2gcJ         gcJ  2gcJ         gcJ 


w'(h+JL     4  *Z  )     4  dw» 
2gcJ         gc7 

-     (w«   4  dw»)    (h  4  dh)   4     (V±dV}2     (wi   *  ^t ) 

2gCJ 

4     -Jt-     (w»   4  dw»)      (Z+d2) 

'•»  ♦  *'       V2       4    gZ       w»  4  dw'h  4  dw»     V2      .     »       q»»z 
2gcJ  «cJ  2-  J   *  *e       T"~ 


fc*cw  ««•  *g<.J       gc 

=     w'h  4  dw»h  4  dh.w»  4  dwdh 

4    J2«l+  2VdVwJ.f  dj£iv£+     2Vdl     dw, 
*«?         2«cJ  2gc'  2gc«> 

4     gjj    (w'^^'Z^'ciZ^wdZ) 

dh.w»   «     2w!VdV     4    jm      wtdZ     =     0 (24) 

2gcJ  «? 
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Momentum  Satiation.     The  net  force  acting  on  the  material 
within  the  differential  control  volume  is  equal  to  the  increase 
of  momentum  flux  of  the  streams  flowing  through  the   control 
volume. 

The  momentum  equation  may  be  written, 

PA  -  (P*dP)A  -  AdZg  {f  *  &£-) 
-  (w*dw)    (V*dV)  -  Vg»  dw-WV 

Vgi  is  the  forward  velocity  component  of  velocity  Vg  with 
which  the  Injected  gas  crosses  the  control  surface.  After 
simplifying, 

-dPA-  AdZ(|*  ♦  &L)     g  -  wdV  ♦  Vdw  -Vgidw 

Neglecting  the  terms  of  second  order 

-AdP-AdZ  g  +  Vg»     dw  .  wdV  -  Vdw (25) 

Definition  of  impulse  function: 

F  .  PA  4  wV 
For  constant  impulse  function  and  constant  ares 

-AdP  .  wdV  *  Vdw 
Momentum  equation  reduces  to 

-A«g  dZ4Vgt     dw.O (26) 

PfTilQBmfgt  fi£  &  la  JCfOss  SLL  *ach  Number. 
T0  s  T(l  ♦  ^M2)        (isentropic  relationship) 

From  the  equation  of  state  and  from  equation  (3) 


^  =  h  s    fi+k*P)RT 


21 


Substituting  the  value  of  T  from  the  stagnation  tempera- 
ture relationship  into  the  isentropic  temperature  relationship. 

U-k  #)RT0  ! 


l*k  M2 


,     Cj.  «    £fU±kJ. (2d) 


'  ^ 


Now 


V  .       AV  .  ?AM  ^/kRT 


s    PAM|     *RTg  v  Substituting  the  value    f> 


1  i+k  h^ 


x         A.M 


kRT, 


1*  1^2 


14k  M2 


where 


C2  =  CX     kRT0 

&  ,     J  <*i*>2     +  a  °T  d(H)2     _     k  d(M) 


»2 


2 

<iw  -     w   *li 
Substituting  the  value  of  w  from  the  above  expression. 

ML 


dw  B 


Lsml     «,     j       2     d(M)2     ,     kdW' 
M2  1*  k;!  M2  i.k  M2 


X         C2  A(M«)*   (1+  ¥**>* 
1+k  Mz 


(29) 


(30) 


(3D 
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Coming  back  to  the  original  values  of  C^  and  C2. 


dw  =     P»(l»k)  .     (1*  11=^2)     «* 

U*kM2)RT0  2  * 


kHTy 
lTljSP 


M2  lfl#  14k  K2 


d(M)2 


k 
14-k 


am[~JL  -        M, 

[3*  4(1+^2) 

■Ma]  d(M)2 ^2) 


DarelooMt  ol  ^UtRWVQfflffg?  Station,  Pflniftfr, 
How  equation   (26)   la  written! 
Vg*   dw  s  APgdZ 

From  aquation  (26),    (27),  and  (32) 

Af     r»U-»fcl  2 dZ 

(1+k  M*)RT0 


-     Vg»        P*U»k) 

(1+k  M2)RT0 


AM   (1+  K|iM2) 


k»t 


UJiclM2 


JL 


2   K2 


♦     X 


i^TmT 


adZ 


M  KRTn 

(1*.  iS^M2^ 


2M2 


Uk  M< 


d(M)' 


Uk  Mz 


jgdZ 


I 


Tk^fo"  .   d(M)2    '  Vg»   :(UWMM 


1+  fcjllf 


k£L 


2M' 


I+2 


14k  M' 


Vg' 


—  (33) 
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In  oraer  to  determine  the  possible  types  of  flow,  it  is 
necessary  to  analyse  the  system  through  consideration  of  the 
second  law  of  thermodynamics.  The  changes  in  entropy  must  be 
positive  because  the  flow  is  (i)  adiabatic  by  assumption  (ii)  ir- 
reversible by  virtue  of  the  mixing  process  of  the  injected  gas 
and  the  ma in- stream  flow. 

Rewriting  equation  (17) 


di  .  k  (1-M2)  d(M)2 

»    2(Uk  M2)  (1  +  ^Jm2) 


From  equation  (17),  with  d§  necessarily  always  positive. 
For  M<1,  d(M)2  must  be  positive. 
For  M>1,  d(K)2  must  be  negative. 

If  the  expression  for  dw  is  examined  by  the  above  criteria, 
it  will  be  observed  that  dw  is  always  positive  for  all  values  of 
the  Mach  Number. 

For  the  analysis  of  expression  (33)  for  vertical  flow,  the 
second  law  of  thermodynamics  imposes  the  following  conditions. 

(i)   d(M)2  is  always  positive  in  subsonic  flow. 

(ii)  d(H)2  it  always  negative  in  supersonic  flow. 

(Hi)  dw  is  positive  for  subsonic  as  well  as  supersonic 
flow. 

After  analysing  expression  (17)  for  the  above  conditions, 
it  is  concluded  that  in  case  of  subsonic  flow  the  only  possible 
flow  is  vertically  upward  flow  with  increasing  Mach  Number. 
Similarly  in  the  case  of  supersonic  flow,  the  only  possible  flow 
is  vertically  upward  flow  with  decreasing  Mach  Number.  The  lim- 
iting value  of  the  Mach  Number  is  unity  in  subsonic  as  well  as 
supersonic  flow. 
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Numerical  Analysis  for  Vertical  Flow  Without  Friction 

The  numerical  analysis  is  done  with  the  help  of  the  momentum 
equation  (26).     In  the  analysis  for  upward  flow,  to  overcome  the 
gravity  effects,  the  injected  gas  must  have  a  forward  component 
of  velocity. 

The  integral  form  of  equation  (26)   is, 

gfdZ  « JVg'   dw 
Assume  Vg*  «  Constant. 
The  Integration  of  the  above  equation  results  in  the  dimension- 
less  parameter,  P*s(Z»  -  Z)    ,  which  is  a  function  of  Mach  Number 

w*  Vg' 

and  was  used  in  the  preparation  of  Table  3,  ?lg*  5,  6. 

ffsyelopaien.t  oj£    f»gU*-Z) 

Vg»w* 

from  equations  (26),    (27),  and    (32) 

g  .    p*u»k)       .    u^m2       .  dZ 

(l*k  M2)RT0  I 


*""'  ?»fe  •  J1^ 


m       4(i+i=lM2)         ^? 


d(M)2 (33f; 


Mm 


*& (»«"'    "-Tie- 

2  ~ 
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And  w*  s    ^*A*V* 

r     ?*A  M*     kRT* 
«     P*   •   A 


\ 


kRT0 


T 


Substituting  equation  (33a)   in  (33b) 

•  A. 


w*     a    p*To 

rHST 


\ 


kRT0 


k+1 


Substituting  equations  (33a),    (33c)   in  aquation  (33*) 


(}*±1) 


g  .     f*R(  2  )(«»!) 
To(l+k  M2)   HT0 


X       2 dZ 


(33b) 


(33c) 


s     AKVg»    «*a ?   ^    2 

AT0     kRT0 


-^fi^f   (i*k)yr*  \jd »  ^m2) 


~ 


♦ k=i. 


4(1+  ^|1m2) 


(1+k  M2)     RTc 


wTTP 


d(M)' 


After  simplification  ana  transposing  the  required  terms  on  the 

left  hand  side  to  get  the  required  form, 

M  


"81      y2(k4l)  d+kfkM2) 


a 


2K 


Hb 


4{U  HM2)  (1+k  M2) 


d(i-0' 


If  the   above  expression  is  integrated  between  M  and  M* 


,/2(k+l)   (u  JsjIm2) 
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^    2  .  /,        Jr. 


4(U*=1M2)  (1*1*2) 


d(M)2 


The  integration  of  the  left  hand  tide  is     ulP+(Z*  -  Z) 

The  author  was  not  able  to  integrate  the  integral  on  the  right 
nana  eide  of  the  equation,  and  a  step-by-step  method  of  numerical 
calculation  was  used. 

Numerical  Analysis  of  Elevation  Effects  Without  Friction 


*1  =  0.5 

M2  s  0.6 

Ti  •  lOOOOg 

?l   -  100  psia.  T0  «  1050°R 

T2  s  0.9523d     T2=  .9523^  x  1050  .  979.48°R  (isentropic  tables) 


ko 


f  1  =  |i-  ■  _iP0 -  0.26995  lbm/ft3 

Ifj    0.37043  x  1000 

Vx  -  49.0172  x  0.5/1000  -  49.0172  x  31.625  x  0.5 

*  1550.1  x  0.5  -  775.03  f .p.s. 

*1  =  ^lvl  s  0.20995  x  775. Od 

■  209.23     lbn/sec. 

T2  3  T0  xO.95238  .  979.4S°R.     (isentropic  tables) 

V2  •  49.0172  x  0.6  x/979.48 

r  31.29  x  49.0172  x  0.6 

s  1533.7  x  0.6  -  920.25  f.p.s. 

|$  r  1.5975,  P*  ■  56.24  psia   (Rayleigh-Line  Tables) 

P2  r  1.5975  x  56.24  «  89.757  psia 
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f  2  «   *9t7S7 »   6?t?6 

1    0.37043  x  979.48     362.825 

«2  «  ? 2V2  »   920.248  x  0,24739 

s   227.66  lbm 
p    s  0.25867  lbm/ft3 
w2-W!  s  18,43  lba 


-  0.24739  lbm/ft3 


Table  3.   Dimensionless  elevation  parometer  for 

subsonic  vertical  flow  without  friction. 


Mach 
Number 

w*  Vg» 

ft. 

z^ 
ft. 

0.5 

123.74xl0"3 

0 

0.6 

74.49xl0"3 

1550.1 

1550.1 

0.7 

39.24xl0"3 

1113.0 

2663.1 

0.8 

12.24x10-3 

724.77 

3387.8 

0.9 

3. 94x10*3 

392.86 

3780.7 

1.0 

0 

120.05 

WOO.  7 

w*  a  251.57   lbn/sec. 

?*  «    0.17351   lbm/ft3 

Vg«  s  700  ft/sec. 

g  9  32.17  ft/sec.2 


to 


Fig, 


5. 


IS  ^pTfi 

Dimension less  ^levacion  parameter  for 
subsonic  vertical   flow  without  friction. r 


-:e*  Lea  -  h) 


's.Mach  Number 


m 


I : 


_ 


_ 


! I ! I I _, E 


Squares  to  the  In<  h 
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Table  4*  Dimensionlesa  elevation  parameter  for  auper- 
aonic  vertical  upward  flow  without  friction, 


f 

Mach 
Number 

w*  Vg« 

^Zft.   : 

Zft. 

; 

: 

i 

2.0 

177.0xlO-3 

0 

0 

1.7 

99.03x10-3 

2455.3 

2455.3 

1.5 

59.10x10-3 

1257.1 

3712.4 

1.4 

40.5x10-3 

590.41 

4302.8 

1.3 

24.3x10-3 

499.69 

4802.5 

1.2 

12.0x10-3 

402.21 

5204.7 

1.1 

2.05x10-3 

271.80 

5475.5 

1.0 

0.0 

99.158 

5574.6 

w*  s  251.57   lbm/sec. 

(?*  -  0.17351   lbm/ft3 

Vg»  ■  700.0  ft/sec. 

g  z  32.17  ft/aec.-2 


JO 


ikjuares  to  the  Inch 
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«c 
Vg»   8  700  f.p.a. 

1  x  0.25387  x  «i  •   700     (18.43) 

32.174 

«  21.756  (13.43) 

=  400.96 

0.25887 
s  1550.1  ft. 

«1  =  0 

i2  r  1550.1  ft. 

Comparison  of  the  Fanno  Line  and  the  Rayleigh  Line  With 
Flow  Having  Constant  Area,  Constant  Impulse  Function, 
and  Constant  Stagnation  Temperature  Fluid  Flow 

In  the  Fanno  Line,  the  stagnation  temperature  and  the  mass 
flow  are  conatant  but  the  impulse  function  decreases.  In  the 
Rayleigh-Line  process  the  impulse  function  and  the  mass  flow  are 
constant  but  the  stagnation  temperature  varies.  In  the  case  of 
the  flow  under  analysis  the  impulse  function  as  well  as  the 
stagnation  temperature  are  constant,  but  the  mass  flow  rate 
varies  through  the  duct.  In  all  three  flows  the  area  of  cross- 
section  of  the  duct  is  constant.  All  these  three  flows  are 
compared  by  plotting  the  following  graphs: 

(1)  P/P*  against  l?/y>*     (Where  V>  is  specific  volume) 

(il)  T/T*  against  S-S* 
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Table  5 

Comparison  of  the  Rayleigh  Line  an 
Line  with  flow  having  Constant   , 
and  Constant  T0, 

d  the  rami 
Constant  F 

0 
» 

Mach 
Number 

■ 
■ 

■ 
« 

:  Fanno 

: 

: 

Line    ; 

* 
* 

: 

Rayleigh  Line  : 

Const.  At  Const. 
F,  and  Const.  T0 

1 : 

|  p/p*  i 
•  „ ,  „;. 

▼A*   : 

• 

P/P* 

1   v/y*  : 
•          * 

P/P*   : 

y/t* 

0.1 

10.94 

0.1094 

2.366 

0.0236 

2.366 

0.5060 

0.2 

5.455 

0.2132 

2.272 

0.0900 

2.272 

0.5240 

0.4 

2.695 

0.4313 

1.960 

0.3137 

1.960 

0.5935 

0.6 

1.763 

0.6343 

1.595 

0.5744 

1.595 

0.7020 

0.8 

1.269 

0.3251 

1.265 

0.3101 

1.265 

0.3410 

1.0 

1.000 

1.0000 

1.000 

1.0000 

1.000 

1.0000 

2.0 

0.4002 

1.6330 

0.3636 

1.454 

0.3636 

1.334 

3.0 

0.2132 

1.9640 

0.1764 

1.533 

0.1764 

2.433 

4.0 

0.1336 

2.1330 

0.1025 

1.641 

0.1025 

2.734 

6.0 

0.0637 

2.2950 

0.0466 

1.630 

0.0466 

3.136 

3.0 

0.0363 

2.3590 

0.0264 

1.6954 

0.0264 

3.234 

10.0 

0.0239 

2.3900 

0.01702 

1.702 

0.01702 

3. 360 

.      1..<»      .  .  . 

2-282 
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Squares  to  the  Inch 
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Table  6.  Comparison  of  the  Rayleigh  Lino  and  the  Fanno 
Lin*  with  flow  having  Constant  A,  Constant  7, 
and  Constant  T0. 


T 


Kach 
Number 


Fanno  Line 


T 


Rayleigh  Line  * 


Const.  A,  Const. 
F,  and  Const.  T0 


T/T* 


_S-2*  |   T/T*  :  .SgS*  |  _SgS* 


T/T* 


0.1 

1.197 

1.363 

0.2 

1.190 

1.085 

0.4 

1.162 

0.4637 

0.6 

1.119 

0.1739 

0.8 

1.063 

0.0392 

1.0 

1.000 

0.0000 

2.0 

0.6666 

0.5247 

3.0 

0.4285 

1.4430 

4.0 

0.2846 

2.3700 

6.0 

0.1463 

3.9740 

8.0 

0.0867 

5.2470 

10.0 

0.0571 

6.2*00 

0.0560 
0.2066 
0.6151 
0.9167 
1.0250 
1.0000 
0.5289 
0.2802 
0.1683 
0.0784 
0.0449 
0.0289 


10.97 
6.325 
1.988 
0.770 
0.1470 
0.0000 
1.2150 
2.7100 
3.9500 
5.8350 
7.2300 
8.3000 


0.2311 
0.2070 
0.1430 
0.0720 
0.0198 
0.0000 

0.4054 
1.2290 
2.1000 
3.6600 
4.9120 
5.9420 


1.197 

1.190 

1.162 

1.119 

1.063 

1.000 

0.6666 

0.4285 

0.2846 

0.1463 

0.0867 

0.0571 


Squares  to  the  Inch 
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Constant  area,  constant  stagnation  temperature,  and  constant 
impulse  function  compressible  fluid  flow  was  investigated  in  this 
report.  In  order  to  achieve  this  type  of  flow  it  was  necessary 
to  have  gas  injected  in  the  main  stream  flow.  In  the  first  sec- 
tion of  the  report  the  ratios  of  pressure,  temperature,  velocity, 
and  mass  rate  flow  through  the  duct  were  derived  as  functions  of 
Mach  Number.  From  these  formulas  a  table  and  a  graph  of  proper- 
ties for  this  flow  were  prepared.  These  ratios  are  Independent 
of  friction  effects  as  well  as  elevation  effects.  In  the  second 
section  of  the  report  the  flow  was  considered  to  be  horizontal 
and  one  in  which  there  was  friction  at  the  wall  of  the  duct.  In 

this  analysis  an  averaging  method  was  used  and  the  results  were 
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expressed  in  the  form  of  the  dimensionless  parameter,    max  . 

D 
This  averaging  process  was  necessary  as  the  author  could  not 
find  an  exact  method  for  the  determination  of  the  flow  proper- 
ties. In  this  analysis  the  gas  injection  had  a  forward  velocity 
component  which  was  equal  to  the  main  stream  gas  velocity.   The 
frictional  coefficient  was  assumed  to  be  constant  throughout  the 
duct.  The  third  section  of  the  analysis  was  made  with  vertical 
flow  without  wall  friction  in  which  the  forward  component  of  in- 
jected gas  was  kept  constant  throughout  the  duct.  In  this  anal- 
ysis the  dimensionless  parameter  a*(Z*  -  Z)  was  developed  and 
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numerical  calculations  were  done  by  a  step-by-step  method,  be- 
cause an  exact  method  could  not  be  found  by  the  author. 

In  all  the  above  cases  it  was  found  that  gas  was  always 


added  to  the  main  stream,  and  that  the  flow  always  tended  to- 
ward a  Mach  Number  equal  to  unity. 

In  the  last  part  of  the  report  a  comparison  of  the  Rayleigh 
Line,  the  Fanno  Line  and  the  flow  having  constant  area,  constant 
stagnation  temperature  and  constant  impulse  function  was  made  by 
plotting  these  processes  on  the  graphs  in  which  coordinates  were 

P/P*  and  v/v*,  and  T/T*  and  S-S*. 
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